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Abstract 

We derive the Lagrangians of the higher-order Painleve equations 
using Jacobi's last multiplier technique. Some of these higher-order 
differential equations display certain remarkable properties like pass- 
ing the Painleve test and satisfy the conditions stated by Juras, (Acta 
Appl. Math. 66 (2001) 25-39), thus allowing for a Lagrangian de- 
scription. 



1 Introduction 

The study of higher-order Painleve equations is interesting from the mathe- 
matical point of view because of the possibility of existence of new transcen- 
dental functions beyond the six Painleve transcendents. In addition such 
higher-order Painleve often have interesting physical and mathematical ap- 
plications. For example it is known that special solutions of equations for 
the Korteweg de Vries hierarchy which are used for describing water waves 
can be expressed via the higher-order Painleve equations. 

The first Painleve hierarchy was first introduced in [Tj. Thereafter many 
results were obtained in the analysis of the higher Painleve equations. Scaling 
similarity solutions of three integrable PDEs namely the Sawada-Kotera, 
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fifth order KdV and Kaup-Kupershmidt equations were considered in [2] 
where it was shown that these fourth-order ordinary differential equations 
(ODEs) may be written as non-autonomous Hamiltonian equations for time 
dependent generalizations of integrable cases of the Henon-Heiles system. 

In j3] it was proved that higher-order members for the first and second 
Painleve hierarchies do not have polynomial first integrals and that their 
solutions can determine new transcendental functions. Lax pairs for some 
equations of these hierarchies are presented in [I] and the Cauchy problem 
for equations of these hierarchies can be solved by an analogy with the Cauchy 
problem of the well known Painleve equations. The Painleve tests for higher- 
order Painleve equations were demonstrated in [5H7]. 

In [8] two new hierarchies of nonlinear ODEs were introduced which were 
called the K\ and K 2 hierarchies and which may be considered as new higher 
Painleve hierarchies. The equations of these hierarchies have all the proper- 
ties that are unique to the famous Painleve equations. 

Shimomura in [9] presented an interesting expression for the first Painleve 
hierarchy which allows us to consider new properties of equations. Poles and 
a - points of the meromorphic solutions of the first Painleve hierarchy was 
studied by Shimomura in [10] , where a lower estimate for the number of poles 
of meromorphic solution is also given. 

In [H] instanton-type solutions and some leading expressions for the sec- 
ond member of the first hierarchy were constructed using multiple-scale anal- 
ysis. Recently Mo in [12] has applied a twistor description of the similarity 
reductions to the case of the KdV hierarchy to obtain the twistor spaces of 
the Painleve I and Painleve II hierarchy. Dai and Zhang [13] have extended 
the results by Boutroux [TH|15] f°r the first Painleve equation to the case of 
the first Painleve hierarchy. The authors have shown that there are solutions 
characterized by divergent asymptotic expansions near infinity in specified 
sectors of the complex plane for higher-order analogue of the first Painleve 
equation. 

Some important results connected with higher-order Painleve equations 
were also obtained in the papers [T6|[T7]. In [16] Claeys and Vanlessen proved 
the universality of the correlation kernel in a double scaling limit near sin- 
gular edge points in random matrix models that were built out of functions 
associated with a special solution of the second member for the first Painleve 
hierarchy. In [T7] the authors established the existence of real solution of 
the fourth-order analogue of the Painleve equation and obtained the solv- 
ability of an associated Riemann - Hilbert problem through the approach of 
a vanishing lemma and found additionally the asymptotics of solutions. 

The Hamiltonian structure of the second Painleve hierarchy was consid- 
ered in [18J. Here the authors introduced new canonical coordinates and 
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obtained the Hamiltonian for evolutions. They also gave an explicit formu- 
lae for these Hamiltonians and demonstrated that these Hamiltonians are 
polynomials in the canonical coordinates. 

The aim of this paper is to obtain the Lagrangians for the four higher 
Painleve hierarchies using the same approach. In recent years much attention 
has been paid to the Lagrangian framework of higher-order differential equa- 
tions. Although a Lagrangian always exists for any second-order ordinary 
differential equation its connection with Jacobi's last multiplier (JLM) [191120] 
is perhaps not very widely known. The credit for resurrecting the JLM, in 
recent years, must go to Leach and Nucci, who have shown how it may be 
used to determine the first integrals and also Lagrangians of a wide variety 
of nonlinear differential equations |21j. While it appears that the connection 
of the Jacobi last multiplier to the existence of Lagrangian functions were 
the subjects of investigation by a few authors in the early 1900's, the precise 
nature of this interrelation was brought out by Rao, in the 1940's [22J. There- 
after it does not appear to have attracted the attention of most researchers 
working in the field of differential equations. 

According to the classical theory of Darboux [23] every scalar second- 
order ordinary differential equation is multiplier variational. The problem 
of finding a Lagrangian for a given ODE is generally referred to as the in- 
verse variational problem of classical mechanics. The necessary and sufficient 
conditions for an equation y" = F(x, y, y') to be derivable from the Euler- 

Lagrange equation |^ — (^f) — 0, was enunciated by Helmholtz [2^] [25] 

in the form of certain identities. 

The variational multiplier problem for higher-order scalar ordinary dif- 
ferential equations has been studied by Fels [26] and Juras [27]. The inverse 
problem for a fourth-order ODE was solved by Fels who investigated scalar 
fourth-order ordinary differential equations of the form 



Fels approach for solving the fourth-order inverse problem was essentially 
based on a modified version of Douglas's [28] classical solution to the multi- 
plier problem as refined by Anderson and Thompson in |29j, who used the 
variational bicomplex theory [30] to derive the multiplier and showed that 
the existence of a multiplier was in a direct correspondence with the existence 
of special cohomology classes arising in the variational bicomplex associated 
with a differential equation. Fels conditions ensure the existence and unique- 
ness of the Lagrangian in the case of a fourth-order equation and it has been 
shown by Nucci and Arthurs [3 1 J and more recently by us [32] that when 



d A u 
dx 1 



f(x,u, 



du d 2 u d 3 u 
dsc doc divC 



)■ 



3 



these conditions are satisfied, a Lagrangian can be derived from the Jacobi 
last multiplier. 

In fact Fels approached the problem using Cartan's equivalence method, 
and arrived at two differential invariants whose vanishing completely char- 
acterizes the existence of a variational multiplier. Unlike the second-order 
case, the multiplier is unique up to a constant multiple. The programme was 
further developed by Juras [27J who studied the inverse problem for sixth and 
eighth-order equations. In fact Juras obtained a similar solution by using, 
however, a more direct approach in the spirit of the variational bicomplex; 
the differential invariants becoming increasingly complicated for higher-order 
systems. By analyzing the structure equations of the horizontal differential 
he uncovered a two-form II with the property dU = mod II, if and only 
if the equation 

d 2n u _ , du d 2n ~ l u, 

dx 2n ' dx' ' dx 2n ~ x 

is multiplier variational. He proved that a Lagrangian, if it exists, is unique 
up to the multiplication by a constant and found functions I±, J 2 , I n , whose 
vanishing provides a necessary condition for the above equation to be varia- 
tional. These functions are not relative contact invariants, but their simul- 
taneous vanishing is a contact invariant condition. 

In [32] the authors made use of the Jacobi Last Multiplier ( JLM) to derive 
Lagrangians for a set of fourth-order ODEs which pass the Painleve test, i.e., 
their solutions are free of movable critical points. Recently the conjugate 
Hamiltonian equations for such fourth-order equations passing the Painleve 
test have also been derived in 1 33 1 - 



2 Four Painleve hierarchies 

Now the first and the second Painleve hierarchy are well known and can be 
written as the following 

N 

^2 t m L m [u} = x, (1) 

m=l 

( d \ M 

( -j- + u J t m L m [u x - u 2 ] - xu - f3 N = 0, (2) 

^ ' m=l 

where N and M are integers, t m , (m = 1, . . . , N) is the sequence of operators 
L m [w] that satisfies the Lenard recursion relation 

d x L m+1 [u] = [dl + 4 u d x + 2 u x ) L m [u] , L [u} = ]-. (3) 
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Taking the operator (J3]) into account we obtain 
Li[u] = u, 



(4) 



L 2 [u] = u xx + 3u 2 , (5) 



L 3 [u] = u xxxx + 10 uu xx + 5 u 2 x + 10 u 3 , (6) 



Z/4[it] u xxxxxx + 14 u u xxxx + 28 u x u xxx + 21 u xx -\- 

(7) 

70 u 2 u xx + 70 u u 2 x + 35 u 4 . 

Using the values of operators L 1; L 2 , L 3 , L 4 and so on we can obtain the 
equations of the first and the second Painleve hierarchies. 

The sixth-order ordinary differential equations of the first and the second 
Painleve hierarchies have the form 

^4 i^xxxxxx 14 u u xxxx -\- 28 u x u xxx -\- 21 u xx ~\~ 70 u u xx -\- 
+70 uu 2 x + 35 u 4 ) + t 3 (u xxxx + 10 u u xx + 5 u 2 x + 10 u 3 ) + (A) 
+ h (u xx + 3 u 2 ) + ti u = x, 

^3 (u xxxxxx 14 u u xxxx 56 u u x u xxx 28 u x u xx 42 u u xx ~\~ 

+70 u A u xx + 140 u 3 u\ - 20 u 7 ) + t 2 (u xxxx - 10 u 2 u xx - (B) 
-10 uu\ + 6 m 5 ) +tx (u xx - 2 m 3 ) - w- /3 3 = 0, 

We see that equations of the first and the second hierarchy have even 
integer orders 2 N — 2 and 2 M respectively. 

Equations (A) and (B) are important and interesting because setting the 
constants t 3 = t 2 = one recovers the Painleve equations. When t± = t 3 = 
these yield equations which we have studied recently. In the case t\ = t 2 — 
they reduce to sixth-order equations which are the third members of the 
first and second Painleve hierarchies. The general case of these equations 
correspond to the first and second Painleve hierarchies. 
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There are two other hierarchies of nonlinear ordinary differential equa- 
tions that have the properties similar to Painleve equations. These hierar- 
chies were introduced in [8] and were referred to in |34J as the K\ and K 2 
hierarchies. These hierarchies can be presented as the following 



N 



(8) 



m=l 



d 



M 



m=l 



U v — - U 



xu — (3 



M 



(9) 



where N and M are integers, t m are parameters of the equation and the 
operator H m may be calculated by means of the formulae 



H n+2 = J[v}Q[v}H n , 
under the conditions 



(10) 



H Q [v} = l, H 1 [v]=v xx + 4:V 2 , 
where the operators Q[v] and J[v] are determined by the relations 



a = D 3 + 2vD + v x , 



D 



dx' 



(12) 



J = D i + 3(vD + Dv) + 2 (D v D~ + D~ v D 2 ) + 



+8(v 2 D- 1 + D^v 2 ), D 



-1 



dx. 



(13) 



Taking conditions f llip and operators fll2p . ( 1131) into account we have the 
operators H 2 and H 3 as the following 



32 

H 2 [v) = v xxxx + I2vv xx + Qv 2 x + — v 3 , 



(14) 
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-^3 M ^uraraxj; "I - 20 V V xxxxxx -\- 60 V xxxxx -\- 134 ^xxxx - ! - 

+136 v 2 r,, + 84 v 2 ^ + 544 v ^ v xxx + 408 v t£„ + 396 v 2 x v xx + 



(15) 



1120 3 cen 2 2 256 5 
u v** + 560 irir; H w 5 . 



Note that hierarchies (jHJ) and (jSJ) can also be presented using another 
operator In terms of this operator these hierarchies take in the form 



N 



J2t k G k [u] = x. (16) 



fc=i 



/ 1 ti \ M 

U- 2 ^C fc [-2 M;E -2u 2 ] -xu-Pm = 0. (17) 

^ x ' fe=i 



Hierarchy f[T6|) can be transformed to (J8J) but hierarchy (TTTj) coincides 
with hierarchy fl9]). The recursion relation G k is determined by the nonlinear 
operator 

G k+2 = J 1 [v]Q[v]G k , (18) 
under the conditions 

G [v] = l, G 1 [v]=v xx + ^v 2 . (19) 

The operator J±[v] takes the form 

J 1 = D 3 + - (D 2 v IT 1 + D~ l v D 2 ) + - (w 2 D" 1 + D~ l v 2 ). (20) 
2 8 

Hierarchies K\ and though similar to the first and the second Painleve 
hierarchies have a fundamental difference in the sense that we cannot trans- 
form equations of hierarchies (fl6|) and (fTTj) to hierarchies ([I]) and (|2J) • More- 
over the hierarchy has even integer order except Gk (k = 1,2,...) and 
hierarchy K 2 also has even integer order except 6 k (A; = 1,2,...). 

The fourth order equation corresponding to the hierarchy K\ takes the 
form 
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h (u xxxx + 12 u u xx + 6 u 2 x + + h ( u xx + 4 u 2 ) = 0. (C) 

At t2 = equation (C) is the first Painleve equation but at t ^ the 
forth order equation differs from the the fourth order equation of the first 
Painleve equation and we hope that this one may give a new transcendal 
function. 

On the other hand the sixth-order equation from hierarchy K-i can be 
written as 

^"2 ^'^ j xxxxxx I ^ x xxxx ^ Vj xxxx I 14 ^-^xx xxx 28 11j 11j x H> xxx 

28 

-28 u 2 x u xx - 21 u u 2 xx - — u u x - 14 u 2 u x u xx + 14 u A u xx + 

4 7 \ & 
+28 u 3 u 2 x - -u 7 j +ti (u xxxx + 5 u x u xx — 5 m 2 

— 5 u u x -\- m 5 ) — X u — 02 = 0. 

Equation (D) is a sixth-order nonlinear ordinary differential equation with 
properties similar to the Painleve equations but cannot be transformed to 
the equation of the second Painleve hierarchy. This equation does not have 
a first integral in the polynomial form and it is possible that it determines a 
new transcendental function. 

In the following section we find the Lagrangians for the nonlinear ordinary 
differential equations (A), (B) and (D). 



3 Inverse problem for sixth-order equations 
and their Lagrangians 

Consider a sixth-order equation in the normal form, u 6 = f(x, u, «i, m 2 , u 3 , m 4 , m 5 ). 
Here we introduce the abridged notation = d k u/dx k . The following theo- 
rem due to Juras gives the necessary and sufficient conditions for a sixth-order 
equation to admit a variational multiplier [27]. 

Theorem. A sixth-order ordinary differential equation admits a varia- 
tional multiplier and non-degenerate third-order Lagrangian if and only if 
following two conditions are satisfied 
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9 du 5 [ x[ du 5 )} *W x[ duJ + 8lW x[ du 5 > 

3 <9-u 5 v <9w 4 3 <9w 5 <9w 4 v <9w 5 ' 3 <9w 3 v <9w 5 ' 3 <9w 5 <9w 3 ' 

_ JL(^) 5 _ Ll^Lf^L _ 1 , 9/ 2 df_ _ 19/9/ _ 9/ 
x ^du 2 243^<9-u 5 ' Tl^du^' du 4 9^du 5 ' du 3 3du 5 du 2 du\ 
and 

U ~3^W 30u B *W X W 

27 <9w 5 ' 3 <9w 5 <9w 4 <9« 3 ' 
Prove. Suppose the sixth-order equation 

u 6 = f(x,u,u 1 ,u 2 ,u^u i ) 

is independent of u^. Then it admits a variational multiplier and a non- 
degenerate third-order Lagrangian if and only if the following two conditions 
are satisfied: 

= DUB ' MB + MB - & and = -2D-4B + ft- 

3.1 The Jacobi Last Multiplier and construction of La- 
grangians for sixth-order equations 

In this section we describe the connection of the Jacobi Last Multiplier with 
the Lagrangian function for sixth-order ODEs. 

Proposition. Let uq = f(x, u, u±, u 2 , w 3 , w 4 , M5) be a sixth-order ordinary 
differential equation which admits a Lagrangian L. Then the function Ai := 

(f^r) , is a Jacobi last multiplier, i.e., it satisfies the equation + J^-M = 
0, where u 5 = u xxxxx . 

Proof : Considering the higher-order Euler operator, E, the Euler- 
Lagrange equation of motion for the ODE w 6 = f(x, u, u±, 115) is given 
by 

EiL ^^- D ^ + D '^- D ^-°- < 21 > 
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where L = L(x,u,Ui,u 2 ,u 3 ) is a third-order Lagrangian. It is obvious 
that the partial derivatives of L, namely L u , L U1 , ....L U3 are all functions of 
x,u, uz- Upon expanding the Euler-Lagrange equation we find that 

= E(L) = U5L U2U3 —[2u5L U3U3X +usL U2U3 +f(x,u,ui,...,U5)L U3U3 

-\-2u^u^L U3U3U3 +2u 5 u 1 L U3U3U +2u 2 u 5 L U3U3U1 +2u 3 u 5 L U3U3U2 ]+teTm8 independent of u 5 . 

Here the subscripts denote partial derivatives with respect to the indicated 
variables. Since the partial derivative of this equation with respect to u§ 
must also be identically zero, we find that 3D X (L U3U3 ) + ^£-(L U3U3 ) = 0. 
Let be M (3) = L U3 u s , then the above equation, E(L) — is expressible as 

be D x ^logM^^ + = 0, showing thereby that the Jacobi Last multiplier 

is given by 

M = Mf 3) . □ 

Remark: Note that for the fourth-order ODE, = f(x,u,...,Us), ad- 
mitting a second-order Lagrangian the analog of ( 13. ip is the following equa- 
tion: 

df 

D x (logM ( 2 2) ) + ^- = 0, 

so that the JLM is A4 = M^L where M( 2 ) = L U2U2 . On the other hand for 
the second-order ODE, u 2 = f{x,u,Ux), it is the solution of 

df 

D x (\ogM {1) ) + ^- = 0, 

with M = M ( i) = L U1U1 . 

Equation (??) provides us a tool for determining the Lagrangian of a 
fourth-order equation once a solution of the defining equation for the JLM 
, Ai, is obtained from (??). In fact in the event / is independent of u^, so 
that the condition (??) is trivially satisfied one obtains the solution Mr$\ = 
constant, which may be set equal to unity, without loss of generality. In such 
a situation the Juras conditions are also considerably simplified as evident 
from the Corollary 1. 

3.2 Determination of the Lagrangians 

We wish to determine a nondegenerate third-order Lagrangian L = L(x, u, Ui, u 2 , u 3 ) 
such that E(L) = 0, where JMj ^ 0, where E is the Euler-Lagrange op- 
erator E = £ - D x (-£-) + £>H(gjy - £>x(gy > and D * denotes the total 
derivative operator D x = + wij^ + u 2 -J^ + ^3^-. If there is a third- 
order Lagrangian satisfying the conditions stated in Theorem 3.1, one says 
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that the ordinary differential equation w 6 = f(x,u,Ui,u 2 ,u 3 ,U4,u 5 ) admits 
a variational multiplier. 

In the new notation equation (A) is given by 

t 3 (u 6 + 14im 4 + 28w 1 w 3 + 21ul + 70u 2 u 2 + 70uuj + 35w 4 ) + 

+t 2 (u 4 + I0uu 2 + hu\ + 10m 3 ) + ti(u 2 + 3u 2 ) = x. 

Proposition. Equation (A) admits a Lagrangian description with La- 
grangian L = t 3 {-\u\ + 7u 5 - 35u 2 ul + 7w«|) + t 2 {\u\ - 5mm? + |u 4 ) + 
ti(— TjM 2 + u 3 ) — x u, where = uS k \ k = 1,2, 

Outline of the proof : In order to show this we will adopt the technique 
used in |32j, to derive Lagrangians for a certain class of fourth-order ODEs, 
namely that of the Jacobi Last Multiplier (JLM). For a sixth-order ODE 
written in the form 

M 6 = ...,U B ) 

one can rewrite the equation as a first-order system: u\ = v, v i = w, w\ = 
s, Si — t, t\ — r, ri = f(x,u,v,w,s,t,r), where the subscript 1 denotes 
differentiation with respect to the independent variable x. Then by definition 
the JLM, A4, for the above system of first-order ODEs is defined as the 
solution of the following equation ^ + (S + f + f; + S + S + f) = 
0. Since in case of eqns (A) and (B) the function / is independent of u$ i.e., 
of r in this notation it follows that dl °^ (3) =0 =^ M( 3 ) = constant. 
Furthermore since the JLM is connected to the Lagrangian by the following 
relation M( 3 ) = where 113 = therefore setting the constant in (??) 
to be —t 3 we have 

d 2 L u 2 

-Q^p = ~h => L = -t 3 y + R(x, u, Ui,u 2 )u 3 + V(x, u, u 1} u 2 ). 

Consequently the determination of the Lagrangian essentially reduces to 
finding appropriate functions R and V such that the Euler-Lagrange equa- 
tion E(L) = 0, reproduces the desired equation, namely (??). Detailed 
calculations however show that it is possible to choose R = so as to sim- 
plify the resulting calculations, and therefore the Lagrangian is of the form, 
La = t 3 (— \u\) + V(x,u,ui,u 2 ). To deduce the unknown function V we 
substitute this form of the Lagrangian into f l2T]) and compare the resulting 
equation with our original sixth-order ODE. One finds that in this case 

dL A „ fdL A \ 2 (dL A \ 3 / dL A 



gives 

-^3^6 = V u — V XUl — u\V UUl — u 2 V UlUl +U4 : V U2U2 +V XXU2 +2uiV XUU2 +2u 2 V XUlU2 +2u 3 V XU2U . 2 
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~l _ '^'l^ / utt«2 ^"l^u\u\ui ~l" U^V U2U2U2 + 2lLiU2 V UUlU2 + 2U\U 3 V UU2U2 4~ 2U2U3 V U \U2U2 ■ 

Inserting the value of uq from the original equation and equating the coeffi- 
cient of W4 we find that 

V U2 U2 = t 3 14u + t 2 

2 

which leads to the following solution, namely V = t 3 7uu\ + t-r^ + Tu 2 + S 
where T and S are functions of x, u, U\. Once again we may set T = so that 
La = ^3 {— \u\ + 7mm 2 ,) + ^2(^2) + S, and it remains therefore to determine 
the unknown function S. From the remaining terms we find that one must 
have S u - S XUl - UiS UUl - u 2 S UlUl = t 3 (70u 2 u 2 + 70mm 2 + 35m 4 ) + t 2 (l0uu 2 + 
hu\ + 10m 3 ) + ti{u2 + 3u 2 ) — x. Next equating the coefficients of U2 it is seen 

2 

that —S = t 3 35u 2 u 2 H-^mm 2 + ti^r + Ku\ + N(x, u). Again choosing K = 0, 
we have ultimately from (13.21) . 

—N u = t 3 35u A + t 2 l0u 3 + h3u 2 - x. 

This yields 

5 

—N = t 3 7u 5 + t 2 -u 4 + t x u 3 - xu, 

and finally the following expression for the unknown function 5* 

1 5 
S = —t 3 35u 2 ul — t25uul + ti(— —u 2 ) + t 3 7u 5 + t2-u 4 + t\u 3 — xu. 



We find finally that the expression for the Lagrangian of eqn. (A) is 
L \ = /.; ( — -m 2 + 7m 5 — 35u 2 m 2 + 7uu\ j +t 2 (1^2 ~~ 5mm 2 + -Uj +ti ( — -ul + u 3 ) — .r //. 



□ 

In a similar manner we find that the sixth-order equations (B) and (D) 
also admit a Lagrangian description, which are stated below. 

Proposition. The Lagrangians associated with the equations (B) and 
(D) are given by 

L B = h [-\ u \ - 2 U j + *2 [\u 2 2 + 5u 2 ul + u 6 \ + 

( 1 5 , \ I > (22) 

+t 3 ( ~2 U l ~ 2 U& ~ 35« 4 «i ~~ 7u 2 ul j - - xu 2 - (3 3 u, 

( 1 7 7 7 7 m 8 \ 

t" / ? ' 9 9 9 4 9 '-i i»t 4 9 \ 

L D = t 2 \ --u 3 + -u x u 2 --u u 2 + ~u x + -u u x - 7u u x - — I + 
+*i (\ u l + l u2 ul ~ \u\ + y \ - X - x u 2 - 2 u. (23) 
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4 Conclusion 



Let us briefly discuss the results of this paper. We have found the La- 
grangians La, Lb and Lp for the sixth-order nonlinear differential equations 
from the Painleve hierarchies. These Lagrangians are generalizations of the 
well known Lagrangians of the Painleve equations. It is interesting to look 
at the properties of the Lagrangians La, L B and L D . For example we know 
that there exists the following symmetry of equation (B) when u(x) — > —u(x) 
and /?3 — > — /?3. This symmetry is preserved for the Lagrangian as well. 
Equations (A), (B) and (D) posses the Painleve property. The Lagrangians 
of several mechanical systems usually involve a difference of their kinetic and 
potential energies respectively. Here this property holds for the Lagrangians 
La and Lb at £2 = £3 = x = 0. However this is not true in the general case 
and as such the Lagrangians derived here are really examples of nonstandard 
ones. The question of irreducibility of the higher Painleve equations is an 
open problem. 
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